Thermometry of ultracold fermions by (super)lattice modulation spectroscopy 
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We theoretically consider non-interacting fermions confined to optical lattices and apply a lattice 
amplitude modulation that we choose to be either homogeneous or of superlattice geometry. We 
study the atom excitation rate to higher Bloch bands which can be measured by adiabatic band map¬ 
ping. We find that the atom excitation rate shows a clear signature of the temperature dependent 
Fermi distribution in the lowest band of the equilibrium lattice as excitations are quasimomentum- 
resolved. Based on typical experimental parameters and incorporating a trapping potential, we find 
that thermometry of one- and two-dimensional systems is within the reach of nowadays experiments. 
Our scheme is valid down to temperatures of a few percent of the hopping amplitude comparable 
to the Neel temperature in interacting systems. 

PACS numbers: 03.75.Ss 05.30.Fk 37.10.Jk 


I. INTRODUCTION 

Ultracold atoms in optical lattices provide a clean re¬ 
alization of various kinds of Hubbard-like models ID¬ 
Es:. These models are well known from the descrip¬ 
tion of solids where they are typically only rough ap¬ 
proximations due to the complex structure of materials. 
Fermionic atoms in optical lattices are of particular inter¬ 
est due to their close analogy to electrons in solids. Ex¬ 
perimentalists have succeeded in observing the fermionic 
Mott-insulator 00 not long after the first Fermi gas 
in a three-dimensional optical lattice was characterized 
0. However, many challenges remain, a major one be¬ 
ing further cooling in order to probe interesting regimes 
such as antiferromagnetic order or unconventional super¬ 
fluidity 0. Moreover, the possibilities to probe ultra¬ 
cold fermionic gases in order to gain precise information 
on their characteristics are still limited. Notably, a reli¬ 
able method to determine the temperature of a fermionic 
gas confined to optical lattices is still lacking. Typically, 
the temperature is measured before and after ramping 
into the optical lattices. For the harmonically trapped 
gas, the temperature can be determined from the inte¬ 
grated density profile obtained in a time-of-flight mea¬ 
surement 0 El- However, this becomes inaccurate at 
low temperatures, i.e. a few percent of the recoil en¬ 
ergy. To determine the temperature within the lattice, 
entropy conservation during the slow loading process into 
the lattice is assumed such that the temperature in the 
lattice can be determined from the initial entropy. This 
is severely limited by non-adiabatic heating processes 
caused by the ramping of the lattice or light scattering 
and fails for in-lattice cooling. Various other schemes 
to directly determine the temperature inside the optical 
lattice have been proposed and some have been experi¬ 
mentally tested. For example, the temperature can be 
determined from the double occupancy which is sensitive 
to thermal fluctuations mm for temperatures on the 
order or above the on-site interaction strength. A tem¬ 
perature measurement has been suggested based on the 
fluctuation-dissipation theorem and a spatially-resolved 


density measurement m which requires in-situ resolu¬ 
tion in the measurement. Other theoretical proposals in¬ 
clude using Raman spectroscopy, transferring the atoms 
to a third hyperfine state, such that the Raman signal 
depends on temperature through the Fermi factor [12| , 
thermometry by light diffraction from atoms in the lat¬ 
tice which results in density-density fluctuations in the 
scattered intensity (T3|, or by studying the response to 
an artificial gauge held [TT] . More recently, the tem¬ 
perature of fermions in an anisotropic three-dimensional 
lattice has been determined from a measurement of the 
nearest-neighbor spin correlator j!5) and by spin-sensitive 
Bragg scattering of light [T5] . All these methods have 
their limitations and most cannot be extended into the 
low-tenrperature regime of interest. Review im gives an 
overview on some of the thermometry schemes. 

In this work, we suggest a scheme how to directly mea¬ 
sure the temperature by means of a time-periodic modu¬ 
lation of the lattice amplitude. Lattice modulation spec¬ 
troscopy is a well established probe of strongly interact¬ 
ing ultracold atoms in optical lattices as first introduced 
for bosonic atoms m- In bosonic systems, energy ab¬ 
sorption imprints a characteristic signal in time-of-flight 
absorption images after sudden switch-off. The energy 
absorbed by the system is typically estimated from the 
broadening of the central peak of the momentum distri¬ 
bution. Thus, precise information on the excitation spec¬ 
trum can be obtained such as a gapped spectrum in the 
strongly interacting regime mm- The same procedure is 
not possible in the case of fermionic systems because here 
the momentum distribution only depends on temperature 
by a smearing of the Fermi surface such that very high 
momentum resolution is necessary in order to extract the 
temperature. However, it was proposed ca that a mea¬ 
sure of the double occupancy gives access to the pairing 
gap or the interaction energy in the Mott insulator for 
attractive or repulsive Fermi gases, respectively, as well 
as information on spin ordering through nearest-neighbor 
correlations. From the theory side, the double occupancy 
as a response to lattice amplitude modulation has been 
extensively studied [201 - E7] . Experimentally, the former 
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was used to identify the fermionic Mott insulator [3] and 
nearest-neighbor correlations have been probed [28]. The 
decay of doublons created by the lattice modulation has 
been measured in order to determine the doublon life¬ 
time [22] • Lattice amplitude modulation has also been 
used as a spectroscopic probe in order to map out higher 
Bloch bands in a quasimomentum -resolved fashion [30] 
or study interband dynamics Ell- 

Motivated by the lack of suitable thermometry schemes 
we propose in this work to use the lattice modulation 
spectroscopy in order to directly determine the tempera¬ 
ture of non-interacting fermions confined to optical (su- 
per)lattices. We consider lattice amplitude modulation 
of homogeneous and superlattice geometry that create 
quasimomentum -resolved excitations that enable us to 
observe signatures of the temperature dependent Fermi 
distribution. While the commensurate modulation con¬ 
serves quasimomentum, the superlattice modulation does 
not which is of interest for many other applications. It in¬ 
jects quasimomentum into the system whenever an exci¬ 
tation is created due to the dimerization of the initial lat¬ 
tice geometry. We show that a temperature determina¬ 
tion is possible in both cases for one- and two-dimensional 
systems and within the reach of current experiments. 
This paper is organized as follows. In Section [il] we in¬ 
troduce the theoretical and experimental framework. We 
define the equilibrium system and the lattice amplitude 
modulation of different geometries. We introduce the de¬ 
tection scheme by adiabatic band mapping and the atom 
excitation rate as our observable in linear response the¬ 
ory. In Section |III| we develop a multiple band model 
of the equilibrium system for typical experimental pa¬ 
rameters and compute the atom excitation rate to higher 
bands by a homogeneous lattice modulation. We inves¬ 
tigate how this serves as a thermometer in presence of 
an external trapping potential. In Section IV we explore 
the possibilities of thermometry by a lattice modulation 
of superlattice geometry applied to both, a homogeneous 
and a dimerized, equilibrium lattice. In all cases, we 
consider the experimental feasibility in detail, such as for 
example requirements on frequency resolution and signal 
strength. 


II. MEASUREMENT SCHEME 



c) superlattice modulation of dimerized lattice 



FIG. 1. The time-dependent amplitude modulation of the 
optical lattice between the two configurations indicated by 
solid and dashed-dotted lines with modulation frequency uo. 
We consider different amplitude modulation schemes. In a) 
the modulation is commensurate with the underlying homoge¬ 
neous equilibrium lattice and conserves quasimomentum. In 
b) the modulation is of dimerized (superlattice) geometry and 
applied to a homogeneous lattice such that quasimomentum is 
not conserved but quasimomentum n/a = kL is injected into 
the system. In c) the superlattice modulation is applied to an 
equilibrium lattice of dimerized geometry such that quasimo¬ 
mentum is conserved. The mean amplitude of the equilibrium 
lattice is indicated by the horizontal dashed line. 


where ^ (x) denotes a fermionic annihilation (cre¬ 
ation) operator and y, is the chemical potential. We 
approximate the optical lattice potentials by Vo (a;) = 
JV = i Vo,i sin 2 (k^Xi) with the laser wave number kL and 
lattice spacing a = Xi denotes the directions x, y 

and z. Tuning the laser intensity in each direction allows 
one to choose different lattice strengths Vo,?: and, conse¬ 
quently, implement different lattice geometries. Three- 
dimensional (3D) cubic lattices, two-dimensional (2D) 
pancakes or one-dinrensional (ID) tubes of atoms can 
be realized. Additionally, we consider an external trap¬ 
ping potential. At the position of the atoms, it can 
be approximated by a harmonic potential, Vr(x) = 
(to/2) to 2 i x 2 . This effective trapping potential in¬ 
cludes the envelope from the lattice beams. 


A. The unperturbed setup 


B. The perturbation 


We study a non-interacting gas of fermionic atoms of 
mass to in two liyperfine states a =f, i an d with particle 
number N = N* + N\ confined to optical lattices. The 
Hamiltonian is given by 


H 0 = J dx dA (x) 

+ V T {x) - \i 4V (x). 


W v * + r„(x) 


In order to create controlled excitations which can be 
used to gain information on the fermionic state, we con¬ 
sider two different perturbations to the optical lattice 
potential. First, we consider standard lattice amplitude 
modulation |]l8j which is commensurate with the under¬ 
lying equilibrium lattice as shown in Fig. [l] a). The 
perturbing potential is given by 5V(x) = sin 2 (fci,x). Sec¬ 
ondly, we consider a superlattice modulation scheme as 
shown in Fig. |T] b). The perturbing potential is then 
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quasimomentum k x 

FIG. 2. The three lowest Bloch bands of a homogeneous opti¬ 
cal lattice along ^-direction and possible excitations by lattice 
amplitude modulation. We consider Vo, x = 7E r in units of 
the atomic recoil energy E r = (hkL) 2 /2m. In the case of 
a commensurate modulation a) the perturbing frequency is 
chosen such that atoms are transferred from the lowest band 
a x = 1 to the first excited band a' x = 2 with conserved quasi¬ 
momentum k x - In the case of a superlattice modulation b) 
the perturbing frequency is chosen such that atoms are trans¬ 
ferred from the lowest band a x — 1 at quasimomentum k x 
to the second excited excited band a' x = 3 at quasimomen¬ 
tum k x + kL which is equivalent to exciting to the shifted 
band E ax= 3 (k x + fei) (dashed line) at same quasimomentum 
k x . The shifted bands E ax= \{kx + kh) and E aj . = 2 (k x + fci) 
are also depicted by dashed lines. In both cases the mini¬ 
mum and maximum value of possible excitation energies (red 
arrows) are related to the initial filling of the lowest band in¬ 
dicated by a thicker line that corresponds to quarter-filling in 
this figure. 


given by SV(x) ~ sin(fcLir). Both perturbations are pe¬ 
riodically modulated with modulation frequency oj and 
small modulation amplitude A. In both cases, we choose 
a modulation of the lattice potential along the direction 
x , regardless of the dimensionality of the underlying equi¬ 
librium system. The full lattice potential becomes time- 
dependent V(x,t) = Vq(x) + Asm(ujt)5V(x). Finally, 
we also apply the superlattice modulation to a dimerized 
equilibrium system as depicted in Fig. |T] c). 


C. The detection scheme 

Determining the energy absorbed by fermionic systems 
from the momentum distribution in time-of-flight images 
is much harder than for bosons. The momentum distri¬ 
bution is step like due to fermionic statistics and not as 
sensitive to heating. In the following, we propose a sim¬ 
ple measurement scheme of the response of the atoms to 
the lattice amplitude modulation. The concept relies on 


the excitation to higher Bloch bands combined with the 
adiabatic band mapping technique. 

We assume the lattice potential to be sufficiently deep 
such that the Bloch bands are well separated and ini¬ 
tially only the lowest band is occupied. We choose lat¬ 
tice modulation frequencies such that we excite from the 
lowest band a x = 1 to the first excited band a x = 2 or to 
the second excited band a' x = 3 along ^-direction. Pos¬ 
sible excitations of the modulation along ^-direction are 
sketched in Fig. [2j In the case of commensurate lattice 
modulation the transition matrix elements are non-zero 
for conserved quasimomentum Afc ss 0 such that excita¬ 
tions occur vertically in quasimomentum space. We con¬ 
centrate on excitations to the first excited band a x = 2. 
We comment on this choice in Sec. IIII Al in more detail. 
For each excitation energy Hco only one initial quasimo¬ 
mentum pair ztk x is on resonance due to the inverted and 
stretched form of the a' x = 2 band. Thus, a quasimomen¬ 
tum resolved excitation is ensured and can be detected 
by a measurement of the transferred band occupation 
(nZ^it)} of the excited band a' x without the necessity 
of quasimomentum resolution within the band. Possible 
excitations strongly depend on the filling of the lowest 
band such that the Fermi occupation which depends on 
temperature will be reflected in the response to the lat¬ 
tice modulation. Moreover, the Fermi dependence of the 
response will be stretched in energy space compared to 
the Fermi distribution in the lowest band such that the 
required frequency resolution needed in experiment to 
resolve the Fermi tail is reduced. Band inversion alone 
leads to a stretching of the Fermi distribution by a factor 
two which is strongly enhanced because the excited band 
is much broader in energy than the lowest band. 

In the case of superlattice modulation spectroscopy ap¬ 
plied to a homogeneous lattice the transition matrix el¬ 
ements are non-zero for a quasimomentum transfer of 
Ak x « kL such that excitations occur vertically in quasi¬ 
momentum space only if we shift the excited band in 
quasimomentum space by A k x = kL- We now concen¬ 
trate on excitations to band a! x = 3 in order to obtain a 
similar band structure geometry as for the commensurate 
modulation scheme as shown in Fig. [2]b). Note, that the 
required experimental frequency resolution is further low¬ 
ered compared to the commensurate modulation because 
of the larger width of the second excited band compared 
to the first excited band. 

In the case of the superlattice modulation spectroscopy 
applied to a dimerized lattice excitations occur with zero 
quasimomentum transfer. Dimerization leads to a fold¬ 
ing of the first Brillouin zone such that the lowest band 
is split into two bands. Excitations occur from the lower 
to the upper band with zero quasimomentum transfer in 
the reduced zone scheme which corresponds to a quasimo¬ 
mentum transfer A k x « kL in the extended zone scheme. 
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D. Linear response theory 


We use linear response theory in order to determine 
the rate of excitations ) for the atoms into higher 

bands a'. We assume the perturbing amplitude A to be 
sufficiently small such that linear response is valid. The 
rate of excitations is given by 

W d,{n ~ ) = wh (i ~£ {IWOMI 2 

1 1 n,m 

x e~ [3En 5{fvjj + E n - E m ) }, (1) 

where Z = exp(—/3I5 n ) is the partition sum, f3 = 

1 /{ksT) is the inverse temperature, |n) is an eigenstate of 
the unperturbed system Hq with corresponding eigenen- 
ergy E n and O = f dx Tj. (x) 5V(x)'d’ IJ (a?) is the per¬ 
turbing part of the Hamiltonian. This formula is known 
as the spectral representation of the response and allows 
for an intuitive interpretation. An excitation is created 
whenever the perturbation O couples many-body state 
| n) to many-body state \m) and conservation of energy, 
Huj = E m — E n , is satisfied. 


III. THERMOMETRY BY COMMENSURATE 
LATTICE MODULATION SPECTROSCOPY 

In this section we investigate the response of the sys¬ 
tem to a commensurate lattice modulation at frequen¬ 
cies corresponding to higher-band excitations. We fo¬ 
cus on excitations to the first excited band a' x = 2 and 
study the response of the homogeneous system as well as 
the trapped system within the local density approxima¬ 
tion (LDA). We show how the Fermi distribution of the 
atoms in the lowest band imprints a clear signal in the 
response such that thermometry is possible for one- and 
two-dimensional equilibrium systems. 


A. The two-band tight-binding model 

In this subsection we derive a tight-binding descrip¬ 
tion of the considered situation. We assume a suffi¬ 
ciently deep optical lattice potential such that the emerg¬ 
ing lowest Bloch bands are well separated. Only the low¬ 
est band a = 1 and one excited band a! x are impor¬ 
tant in our considerations and thus taken into account. 
We consider an atom trapped in the d-dimensional pe¬ 
riodic potential Vq(x ) in the absence of an additional 
external trapping potential. We compute the eigenfunc¬ 
tions which are the Bloch functions along the different 
directions <j)aj. (x,) with corresponding eigenvalues (the 
Bloch bands) E aj . (k Xi ) where a Xi denotes the band in¬ 
dex and the quasimomentum k Xj lies within the first Bril- 
louin zone ] — fez,, Ay]. The full spectrum is given by 
E a (k) = Y^i=i E a x (k Xi ) as we consider non-interacting 





quasimomentum k x 


FIG. 3. a) The transition probability |M Q! ^ = 2 (fe a; , fc(,)| 2 for 
exciting atoms at quasimomentum k x from the lowest band 
a = 1 to the first excited band a' = 2 for all k x by com¬ 
mensurate lattice modulation, b) The transition probability 
| M a i (k x )| 2 for exciting atoms to the first a 1 = 2 or second 
a' = 3 excited band at quasimomentum k x = k x by com¬ 
mensurate lattice modulation. We consider a lattice depth of 
Vo, x = 7E r . 


fermions for which the potential is separable. The in¬ 
dex a = {a x ,a y ,az} labels the band. Note, that the 

Bloch functions factorize <^(x) = nti (xi) in the 
tight-binding approximation neglecting the coupling to 
other Bloch bands and the overlap of different sites. The 
unperturbed Hamiltonian (for Vt = 0) becomes 

#o ~ ^2 \ E <x(k) - A 4 

Oi,k,cr 



where c^'i are the fermionic annihilation (creation) op- 
erators. 

In order to construct the perturbing Hamiltonian, we 
need to determine the transition matrix elements of the 
perturbing potential SV(x) in the chosen Bloch basis. 
Initially, only the lowest band a = 1 = {a x = 1, a y = 
l,a z = 1} is occupied. Thus, we are only interested in 
the transition matrix elements given by 


M, 


k'* , 


(a=i,£)->(«',*) = cT dx {x)8V{x)(j) a l = i(as) 


X . ,oc Xi =1 (3) 

*# 1 


where Q x = (L — l)a is the system size in x-direction with 
L the number of lattice sites and x m j n = —(L/2 + l)a 
and x max = La/2. We used that the Bloch functions are 
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orthogonal, i.e. (l/Cl Xi ) dxi (xi) = 

5 a > ,« x . .,k x - We insert the perturbing potential 
5V(x) = sin 2 (kL,x) and find that the perturbation, to a 
good approximation, only couples momenta k and k' ~ k 
as shown in Fig. [3] a). All other matrix elements for 
A k ^ 0 are strongly suppressed. This justifies the ap¬ 
proximation that the quasimomentum is conserved by 
the commensurate lattice modulation. We approximate 


M. 


(aL=l,k)—>(a , ,k') 


— Ma' x (k x , k x ) ,fe Kj 

& 1 


where excitations occur to the band a' = {a' x ,a y = 

f; a z 1} * 

In the tight-binding approximation the perturbing part 
of the Hamiltonian becomes 


Hpert ~ A sin(wt)0, 

O = ^ M a <_ (. k x ) ( c ^,£ CT c a= ifc CT + h.c.^ (4) 

k,a 

The probability | M a i (k x )| 2 of exciting atoms from the 
lowest band to the band cJ x = 2 or a x = 3 with A k x = 0 
is shown in Fig. [3] b). Excitations to the first excited 
band are prohibited at the center k x = 0 and the border 
k x = fcz, of the Brillouin zone due to the symmetry of 
the Bloch functions. The maximum probability lies in 
between. This is advantageous as this region coincides 
with the region of interest around the Fermi surface of a 
half-filled lowest band such that the temperature depen¬ 
dence of the Fermi surface can be probed by the response. 
Moreover, the amplitude is sufficiently strong such that 
enough atoms for a detectable signal are excited at small 
perturbing amplitudes and reasonable perturbing times. 
We supply numbers for an example in Section |III C| in 
presence of an external trap. In contrast, excitations 
along x to the second excited band a' = 3 are non-zero 
for all momenta and the amplitude is about three times 
larger. The transition probability \M 3 (k x )\ 2 is maximum 
at k x = 0 and k x = k^. The corresponding contributions 
dominate the response at the corresponding frequencies 
if considering the full range of excitations. This makes it 
more favorable to use excitations to a' = 2 for a temper¬ 
ature measurement at intermediate lattice heights. Note 
that stronger transition amplitudes to a' = 3 than to 
a' = 2 are in agreement with the expectations for a 
large lattice height where lattice wells almost decouple. 
Approximating each well by a harmonic oscillator, the 
modulation corresponds to a frequency modulation of 
the quantum harmonic oscillator, ojho = x>ho( 1 + p), 
where p is a small parameter. This modulation couples 
to the second excited band, but transitions to the first ex¬ 
cited band are prohibited by symmetry. We confirm that 
\M 2 (k x )\ 2 decreases for increasing lattice depths which 
corresponds to approaching the harmonic oscillator limit 
in a single well. 



FIG. 4. Excitation rate of the atom number to the first ex¬ 
cited band a x = 2 by the commensurate lattice modulation 
spectroscopy of the homogeneous system for a) half-filling and 
b) quarter-filling for temperatures ksT = 0, 0.5J, 1J and 1.5J. 


B. The response of the homogeneous system in one 
dimension 

We start our discussion of results by the simplest sys¬ 
tem which consists of homogeneous ID tubes oriented 
along x-direction. Assuming that only a x = 1 is ini¬ 
tially occupied, it can be shown that the response func¬ 
tion Eq. |T]) in the continuum limit becomes 

1 ' t t t r k L 

— dt{ nZ^ D = jr^J Q d kx |M a , (Ml 2 

x/ [E aiD=1 {k x ) - p)\ S {hu - [E a > x {k x ) - E ax=1 (k x )\ } . 

This expression strongly depends on temperature 
through the Fermi function f(E — /i) = 1/{1 + 

exp [(E — p)/ksT]}. We show the results for excitations 
to the first excited band a' x = 2 for different fillings and 
temperatures in Fig. [4j Here we consider 40 A'-atoms 
and a laser wave length of A = 1064nm. All energies 
are conveniently expressed in units of the atomic re¬ 
coil energy E r = h 2 /2m\ 2 . We use these parameters 
throughout the remainder of this work. Temperatures 
we express in units of the tunneling amplitude of atoms 
in the lowest band between adjacent sites of the optical 
lattice J = (£a x = l.max ~ ^a x = l,min)/4 « 0.04E;r. The 
zero-temperature response vanishes for excitation ener¬ 
gies below the minimum possible excitation energy. At 
zero temperature, the minimum excitation energy corre¬ 
sponds to the excitation energy of an atom located at 
the Fermi surface (cf. Fig. [2] a) ). Smaller energies cor¬ 
respond to states located above the Fermi surface which 
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are empty. At higher excitation energies the response 
is non-zero because atoms below the Fermi surface can 
be addressed. The maximum possible excitation energy 
is always given by the energy needed to excite an atom 
located at k x = 0. At finite temperatures atoms get ther¬ 
mally excited around the Fermi surface. The Fermi step 
softens according to the Fermi function such that the ex¬ 
citations depend on the Fermi distribution in the lowest 
band. The minimum excitation energy is no longer sharp 
but the Fermi dependence is reflected in the response as 
a function of excitation energy. We compare half-filling 
in Fig. [4^) and quarter-filling in Fig. [4jo) for tempera¬ 
tures ksT = 0, 0.5J, J and 1.5J. For both fillings the 
temperature dependence of the response is clearly visible. 
At higher temperature the Fermi step smears out more 
which results in a broadening of the low frequency tail, 
in particular at low fillings. The possible broadening is 
limited by the upper band edge, limiting the resolution at 
high temperatures. The position of the Fermi tail shifts 
to smaller energies for larger fillings due to the depen¬ 
dence of the excitation energy Huj on the quasimomentum 
k x . Furthermore, the Fermi tail is broadened in energy 
due to the enlarged bandwidth of the first excited band 
with respect to the lowest band. The different response 
functions can be clearly distinguished for the displayed 
temperatures. 

C. The response of the trapped system in one 
dimension 

As we have seen that the response depends strongly 
on the filling, we continue our discussion by adding an 
external harmonic trapping potential as present in many 
experimental setups. We treat the trap within the LDA. 
The space-dependence of the trap Vt(x) = V t (x/a) 2 is 
incorporated into the chemical potential. The chemical 
potential thus varies when moving through the trap, i.e. 
n(x) = hq — Vr(x), where ^0 is the chemical potential 
at the center of the trap. We obtain the mean response 

dt(nZ*) 1D /\A\ 2 of the trapped system by summing the 

homogeneous response per lattice site d t (nZ x ) 1D / (|A| 2 L) 
over all chemical potential present in the trap with the 
correct weight and normalizing by the total number of 
atoms. The mean response of the trapped system be¬ 
comes 

d t (r&) 1D 1 [2J p d t (n^y D 

|AP NL\A\ 2 V V* loo ^ ’ 

with /2 = n/2J and V t = (m/2)wf x a 2 . The chemi¬ 
cal potential in the trap center /iq is determined by the 
total atom number N. The mean response for an ini¬ 
tial atom number N = 60 and a trapping frequency 
uj t .,x = 27r x 24.5Hz= 5.6 x 10 ~ 3 E r in a ID system is 
shown in Fig. [5] a) for temperatures ksT = 0.5J, 1J and 
1.5J. Error bars correspond to 5% uncertainty on the ini¬ 
tial atom number N. The results resemble the response of 



hw [E r ] 



FIG. 5. Excitation rate of the atom number to the first ex¬ 
cited band a x = 2 of the commensurate lattice modulation 
spectroscopy of the trapped system. In a) the underlying equi¬ 
librium lattice is one-dimensional with atom number N = 60. 
In b) the underlying equilibrium lattice is two-dimensional 
with atom number N = 4000. In both cases, the response 
is shown for temperatures ksT = 0.5J, 1J and 1.5J with er¬ 
ror bars corresponding to 5% uncertainty on the initial atom 
number N. The vertical lines indicate huj = 3.7 E r for which 
we determine the number of excited atoms (see main text). 
The inset in b) shows the first Brillouin zone of the square 
lattice system. At zero temperature and half filling all states 
below the Fermi surface are occupied (gray) and all states 
above the Fermi surface are empty (white). At a fixed k x 
(vertical black line) different k y contribute such that the re¬ 
sponse is a superposition of constant contributions in the bulk 
and temperature dependence mainly comes into play near the 
Fermi surface. 


the homogeneous system and a temperature dependence 
is clearly visible. However, the Fermi tail is less distinct 
as we sum over different fillings present in the trap. In 
order to resolve the Fermi tail and distinguish the curves 
of different temperatures we estimate from the horizontal 
distance of the curves in Fig. [5]a) that a frequency resolu¬ 
tion A (Hui) ss ±0.02511). is required in experiment which 
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corresponds to a perturbing time t = 27r/Aw ss 9ms. We 

determine the total number of atoms ( nZ x (t )) excited 
during the time t which is the observable in experiments. 
The number of excited atoms in d-dimensional systems 
is given by 

{nZ x (t)} = N x d t (nZ’ B ) dD x t. (6) 

We choose a sufficiently small amplitude modulation of 
A = 0.05V o , x where V$ jX = 7 E r and consider the point 
hu> = 3.7 E r as an example, indicated by a vertical line 
in Fig. [5] a). Using Eq. ([6]) we determine the num¬ 
ber of excited atoms (nS‘ I_2 (t)) to lie between 3 and 7 
for temperatures ksT = 0.5J, 1J and 1.5J. Considering 
~ 100 parallel ID tubes of roughly equal filling in exper¬ 
iment, this gives a difference in number of excited atoms 
of ~ 150 between the chosen temperatures which is mea¬ 
surable by current experimental means. 


D. The response of the trapped system in higher 
dimensions 


We may also probe higher-dimensional equilibrium sys¬ 
tems by lattice amplitude modulation along one direc¬ 
tion. The response of d-dimensional lattices becomes 

d /» 

/ d d k \M a ,(k x )\ 2 


dt(nZ) dD = 

\A \ 2 "> \A*L/ JBZ 

xf Ei(k) — n) 5 {ftw - [E a ^(k x ) - E ax=1 (k x )]} 


L 

2kL 


where the integral runs over the first Brillouin zone (BZ). 
Due to the modulation along the ^-direction only, the res¬ 
onance condition in this expression is set by the change of 
the energy in x-direction, E a t ( k x ) — E aa .-i(k x ), and the 
matrix element depends only on k x . The difference to the 
system of ID tubes enters in the initial Fermi distribution 
where the energy of the d-dinrensional system E\(k ) oc¬ 
curs. For d > 1 the signal remains k x resolved, however, 
many different fcy-quasimomenta contribute at each fixed 
k x . These points have different location in fc-space with 
respect to the Fermi surface such that we detect a super¬ 
position of different points in the Fermi distribution. We 
illustrate this for half-filling in a 2D system in the inset 
of Fig. mb). Consider an intermediate fc x -value as indi¬ 
cated by the vertical line. Along this line a temperature 
dependence will only show up close to the Fermi surface 
whereas the bulk dominantly contributes with a constant 
value to the response. Consequently, the effect of tem¬ 
perature in the response is less pronounced. However, 
in 2D a temperature measurement is still possible but 
requires better frequency resolution in experiment than 
the ID case. Within LDA we obtain 


dt{nZ') 2D 


\A\ 2 


NL 2 \A/- 


4 J 

Vt 


rv 0 


= ) / d,i3 f (n“') 2D (,i)(7) 


where /i = /x/4 J and V t = (m/2)^Ju/ x ijj/ y a 2 . The result 
for trapping frequencies ut,x = 27 tx 24 . 5 Hz= 5.6xlO _3 -E r 


and uit, y = 27r x 29.7Hz= 6.7 x 10 i E r and an ini¬ 
tial atom number N = 4000 are shown in Fig. m b) 
for temperatures k^T = 0.5J, 1J and 1.5J. Error bars 
correspond to 5% uncertainty on the initial atom num¬ 
ber. In order to resolve the Fermi tail and distinguish 
the curves of different temperature a frequency resolu¬ 
tion A (Hco) ps iO-OlE/ is required in experiment which 
corresponds to a perturbing time t ps 23ms for the same 
system parameters as in Section III C We choose a small 
amplitude modulation of A = 0.02Vo and consider the 
point Hu = 3.7E r as an example, indicated by a verti¬ 
cal line in Fig. [ 5 ] b). Using Eq. ([6]) we determine the 
number of excited atoms (n£; =2 (f)) to lie between 150 
and 180 for temperatures kgT = 0.5J, 1J and 1.5J. This 
gives a difference in number of excited atoms of ~ 10 
between the chosen temperatures. This may be demand¬ 
ing to measure but achievable with current experimental 
techniques. 

Note that applying the perturbation along several di¬ 
rections does not help in order to gain on resolution. 
The modulation decouples into the different directions, 
for example a lattice modulation along two directions 
SV(x) = 5V(x) + SV(y ), such that it has the effect of 
two ID perturbations along x and y , respectively. As 
one has the contributions of both directions, the total 
number of excited atoms is only increased by a factor of 
two for a completely isotropic setup. We do not benefit 
from this factor as we are limited in minimum time by 
frequency resolution. 


IV. THERMOMETRY BY SUPERLATTICE 
MODULATION SPECTROSCOPY 

In this section we study the possibilities of thermom¬ 
etry by superlattice modulation spectroscopy. We con¬ 
sider two different equilibrium situations. First, we con¬ 
sider the homogeneous optical lattice system, where all 
wells are equal and of spacing a = A/2, subjected to a 
superlattice modulation. In this situation the symmetry 
of the perturbation is different with respect to the equi¬ 
librium system. The periodicity in space is doubled and 
now 2a whereas the Brillouin zone has half the size with 
respect to the equilibrium lattice. Secondly, we consider 
a dimerized equilibrium lattice such that the superlattice 
modulation scheme has the same symmetry as the under¬ 
lying lattice geometry. 

In the first case, we discuss higher-band excitations along 
the lines of the previous Section |TTT| where we now choose 
to consider the excitations to the second excited band 
a' = 3. The lattice geometry becomes dimerized which 
strongly affects the nature of excitations, namely excita¬ 
tions with quasimomentum k l are created. 

In the second case, we probe a dimerized equilibrium sys¬ 
tem (superlattice geometry) such that the superlattice 
modulation is commensurate with the initial geometry 
[cf. Fig. 0 c)] and conserves quasimomentum. In this 
case, dimerization leads to the emergence of two gapped 
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quasimomentum k x 


FIG. 6. a) The transition probability | M^ a=1 £/)| 2 

for exciting atoms to the second excited band a' = 3 for all 
k' x by superlattice modulation, b) The transition probability 
\M a , ( k x )\ 2 for exciting atoms from the lowest band a = 1 
at quasimomentum k x to the lowest band a' = 1 or second 
excited band a! = 3 at quasimomentum k x = k x + kL- 


energy bands such that excitations occur vertically in 
quasimomentum space from the lower to the upper band. 
As for the lattice modulation spectroscopy we develop a 
tight-binding description for both situations. 


A. The response of a homogeneous optical lattice 

We start our discussion by the first case, considering 
an equilibrium system which is given by a homogeneous 
optical lattice setup and apply a superlattice modula¬ 
tion. We consider the multiple-band model developed 
in Section fill A| As before, the equilibrium Hamiltonian 
H 0 is given by Eq. ([2]). As we consider the superlattice 
perturbation, we need to compute the transition matrix 
elements in Bloch basis given by Eq. <© where we now 
insert the perturbing potential 6V(x ) « sin (kLx) which 
is incommensurate with the underlying equilibrium lat¬ 
tice. We find, that to a good approximation [cf. Fig. [6] 
a) for a' = 3], the superlattice perturbation only couples 
quasimomentum k x in the lowest band to quasimomen¬ 
tum k x + Jcl in higher bands, 

i^l 

~ M a ' x {k x )5p ,fc+(fc^,0,0)• 




FIG. 7. a) A superlattice perturbation can lead to excitations 
within the lowest band E(k x ) with quasimomentum transfer 
A k x = k l. Possible excitations depend on the filling indi¬ 
cated by the chemical potential /r (dotted line), b) The cor¬ 
responding temperature-dependent response of the atoms to 
a superlattice amplitude modulation for a quarter-filled band 
and different temperatures. 


The perturbing part of the Hamiltonian becomes 
ffpert = Asin(wf)0, 

O = M a / ( k x ) (o' - „ c + h.c.^) . 

Z_^ x > \ a'k+(k L ,0,0)(T a=lka J 

k,cr 


In contrast to the commensurate modulation [cf. Eq. 
©]. transitions within the lowest band are possible due 
to the finite quasimomentum transfer, provided that 
empty states at k x + kj, are available. The probabil¬ 
ity | M a ^(k x )\ 2 of exciting atoms within the lowest band 
a = 1 to a' = 1 or from the lowest band a = 1 to the 
second excited band a! = 3 is shown in Fig. |b]b). For 
the probability of excitations within the lowest band we 
find \Mi(k x )\ 2 oc sin 2 {k x a). The excitation probability to 
the second excited band \M 3 (k x )\ 2 is similar in symmetry 
to \Mi{k x )\ 2 but not a pure sinus due to the difference 
in shape of the excited band. In both cases, prohibited 
transfers are located at k x = 0 and k x = kL as for the 
commensurate modulation when exciting to the first ex¬ 
cited band. The response becomes 


dt(nZ) dD = tt (L_ 

\A\ 2 h \2k L 

xf E a -i(k) — n 5 [hio — 


d d k | M a , x (k x )\ 2 

iZ 

\Ea r x {k x T k L ) E ax —i 


( 8 ) 
(M]} > 


in d = 1, 2 or 3 dimensions. 

We now discuss the purely one-dimensional model and 
solely excitations within the lowest band (a' = 1). Atoms 
with quasimomentum k x are excited within the lowest 
band to an unoccupied state k x + if available while the 
system absorbs the energy hu = E(k x + — E(k x ) = 

—2 E(k x ). We approximate E{k x ) = —2J cos(k x a) where 
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FIG. 8. The response to a superlattice modulation along one 
direction of the trapped system within LDA. In a) the under¬ 
lying equilibrium lattice is one dimensional with atom number 
N = 60. In b) the underlying equilibrium lattice is two di¬ 
mensional with atom number N = 4000. In both cases, the 
response is shown for temperatures fcsT = 0.5J, 1J and 1.5J 
with error bars corresponding to 5% uncertainty on the initial 
atom number N. 


J is the tunneling matrix element of the lowest band in 
^-direction. Possible excitations may be represented by 
transferring atoms at constant quasimomentum k x to the 
inverted band E(k x +k B ) as shown in Fig. [7]a). The atom 
excitation rate [Eq. ([sf] can be reduced to 

frdtinZ*- 1 ) _ max(|Mi(fc x )| 2 ) f / Huj\ 

L\A\i 4J V \Jj) 



The temperature-dependence of the response function is 
shown in Fig. [7]b). The Fermi tail in the response is 
stretched by a factor of two with respect to the initial 
Fermi occupation of the lowest band due to the band in¬ 
version. This effect is strongly enhanced by higher-band 


excitations that hence are in the focus of our discussion. 
Consequently, in the following we restrict the presenta¬ 
tion of results for the atom excitation rate to the second 
excited band a' = 3 in the presence of an additional ex¬ 
ternal trapping potential. We study the response to su¬ 
perlattice modulation of the trapped system in the LDA 
along the lines of Sections |III C| and |III D| choosing the 
same system parameters and the same initial atom num¬ 
ber. The mean atom excitation rate in the trap is given 
by Eq. @ in ID and Eq. 0 in 2D where we insert the 
response to the superlattice modulation of the homoge¬ 
neous system [Eq. ([8]) in ID and 2D, respectively]. The 
total atom excitation rate to the second excited band 
of the ID and 2D system with 5% uncertainty on the 
initial atom number is shown in Figs. [8] a) and b), re¬ 
spectively. We observe better frequency resolution com¬ 
pared to the case of commensurate modulation due to 
the enhanced width of the second excited band. In or¬ 
der to resolve the Fermi tail and distinguish the curves 
of different temperature in ID a broad frequency reso¬ 
lution A (hui) ~ ±0.15-EV. is already sufficient in experi¬ 
ment which corresponds to a perturbing time t ss 2.5ms 
for the chosen parameters. We consider a typical time of 
t = 10ms and a perturbation amplitude A = 0.05Vo,®- At 
the point Huj = 7E r indicated by a dashed line in Fig. |8]a) 
as an example. Using Eq. © we determine the number 
of excited atoms (n% =3 (f)) during time t to lie between 
2 and 6 for temperatures k B T = 0.5J, 1J and 1.5J which 
shows that we excite a couple of hundreds of atoms when 
considering ~ 100 parallel ID tubes of approximately 
equal filling with an atom number difference of ~ 150 
which is detectable in experiment. In 2D the improved 
frequency resolution is a clear advantage compared to 
the commensurate lattice modulation. We find that a 
frequency resolution A (Huj) ss ±0.025E? r is already suf¬ 
ficient in experiment which corresponds to a perturbing 
time t ~ 9ms. Considering an amplitude A = 0.05Vo,a; 
at Huj = 7E r we excite atom numbers {n%^ 3 (t)) be¬ 
tween 300 and 400 for temperatures k B T = 0.5J, 1J and 
1.5J with an atom number difference of approximately 30 
atoms which is easier to detect in experiment than the 
smaller atom number difference in the case of commen¬ 
surate lattice modulation. 


B. The response of a dimerized lattice 

In this subsection we consider a one-dimensional 
dimerized system as the equilibrium system and apply 
the superlattice modulation which has the same geome¬ 
try as the equilibrium system. The equilibrium superlat¬ 
tice is engineered to have a constant bottom offset but 
alternating lattice height [cf. Fig. |T] c)]. The equilibrium 
tight-binding Hamiltonian describing the considered su- 
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4 

huj[J\ 


FIG. 9. a) The superlattice modulation of the dimerized (su¬ 
perlattice) equilibrium lattice ( J' = 2 J) leads to excitations 
from the lower energy band — Ek to the upper energy band 
+Ek at conserved quasimomentum k. Possible excitations de¬ 
pend on the filling given by the chemical potential p (dotted 
line), b) The corresponding temperature-dependent response 
function to a superlattice amplitude modulation for a quarter- 
filled band and different temperatures. 


perlattice configuration is given by 

H q k-J ( c ],* c i+i,<r + h - c 0 

j odd,cr 

— J 'y ' ( C i.g-Cj+ l,g ~t~ h.C.) — fl n j,cri (9) 

j even,(j j,a 

where hopping amplitudes J' > J > 0 alternate between 
neighboring sites. The unit cell contains two sites that 
may be assigned to two sub-lattices. Consequently, the 
first Brillouin zone has half the size ] — fcz,/2, &l/ 2] com¬ 
pared to the homogeneous lattice. We employ a Fourier 
transform combined with a Bogoliubov transformation 
(cf. Appendix |A| for details) and obtain 

Hq = Ek{a} k G a.k,<j — Pk tCr Pk,cr) 
k,a 

- V'52( a k,cr a k,cr + Pt t<r Pk,<r), ( 10 ) 

k,a 

which leads to the emergence of two energy bands ±E)~ = 
± \J J’ 2 + J 2 + 2 JJ’ cos(2 ka) gapped by an energy A = 
2(J' — J) and of total band width W = 2 (J' + J) [cf. 
Fig. [9] a)]. Here k = k x denotes the quasimomentunr in 
x-direction. The operators ct^\ and /3^ annihilate (cre¬ 
ate) quasiparticles in the upper or lower band, respec¬ 
tively. Excitations due to the superlattice modulation 
occur from the lower band — Ek to the upper band +Ej. 
at conserved quasimomentunr k as shown in Fig. EM- 
This is in contrast to the case of a homogeneous equi¬ 
librium lattice where the superlattice modulation intro¬ 
duces quasimomentum A k = kL into the system. How¬ 


ever, here, quasimomentum is conserved as the equilib¬ 
rium lattice is already dimerized and A k = k^ = 0 in 
the reduced zone scheme of the superlattice. For the 
atom excitation rate to the upper band we obtain (see 
Appendix |A| for details) 

hd t (nl PP erband ) _ \C(w)\ 2 F{u>,P,n) 

L|A | 2 1 j 


where the functions |C(w)| 2 and \g'(oo)\ do not depend 
on temperature (cf. Appendix [A]) and the temperature- 
dependent part is given by 


F(u,/3,n) = 


e /3u/2 _ e -/3w/2 


e~^(e^ + e^“/ 2 )(e^ + e~P^/ 2 ) 


The strong temperature-dependence of the response is 
shown in Fig. |9]b). The temperature-dependent tail is 
also stretched by a factor two in energy as for the ho¬ 
mogeneous equilibrium lattice when exciting within the 
lowest band but excited states are now clearly separated 
from the initial distribution which is an advantage. Con¬ 
sequently, a simple thermometry scheme is available for 
the dimerized equilibrium lattice by superlattice modu¬ 
lation spectroscopy by counting the number of excited 
atoms in the upper band if the parameters allow for suf¬ 
ficient frequency resolution. Here, we do not consider 
explicit experimental parameters as our aim was to out¬ 
line the conceptual scheme in in the dimerized superlat¬ 
tice. However, it can be extended, incorporating typical 
experimental parameters. 


V. CONCLUSION 

In this work we have studied a non-interacting Fermi 
gas confined to optical lattices. We probed the response 
to a time-periodic modulation of the lattice amplitude. 
We investigated the temperature-dependent atom excita¬ 
tion rate to higher Bloch bands and demonstrated that 
it shows clear signatures of the Fermi distribution of 
the equilibrium system. We explored the possibilities 
of thermometry for different geometries of the equilib¬ 
rium system as well as for different amplitude modula¬ 
tion schemes, either commensurate or incommensurate 
with the equilibrium system. 

We first considered a homogeneous equilibrium lattice 
subjected to both a homogeneous lattice amplitude mod¬ 
ulation [Fig. [I] a)] and a lattice amplitude modulation of 
superlattice geometry [Fig. l] b)]. We studied the exci¬ 
tations to higher Bloch bands on the basis of a multiple- 
band model which we constructed for typical experimen¬ 
tal parameters. We found that for the homogeneous mod¬ 
ulation quasimomentunr is conserved and excitations to 
the first excited band are suitable for thermometry. In 
contrast, the superlattice modulation transfers quasimo¬ 
mentum kL into the system and we consider excitations 
to the second excited band. In both cases, the response 
displays a clear signature of the temperature-dependent 
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Fermi distribution in the lowest band for one -and two- 
dimensional equilibrium systems where we considered 
temperatures on the order to a few percent of the hopping 
amplitude J. The Fermi distribution is strongly broad¬ 
ened in energy due to the much larger bandwidth of the 
higher bands compared to the lowest band. We benefit 
from this because the experimental frequency resolution 
required to resolve the Fermi dependence is achievable 
within typical durations of the perturbation. We estab¬ 
lished close relation to experiments by incorporating an 
external trapping potential and by estimating the num¬ 
ber of excited atoms in order to verify the measurability. 
Additionally, we investigated the response of a dimerized 
equilibrium system subjected to the superlattice modu¬ 
lation setup [Fig. [l] c)]. We found that this setup is also 
suitable for thermometry due to the emergent two-band 
structure. 

We emphasize that the temperature-dependence of the 
atom excitation rate becomes more distinct for decreas¬ 
ing temperatures such that our scheme extends down 
to temperatures that have not been reached in experi¬ 
ment so far and covers the regime of interest below the 
Neel temperature where antiferromagnetic ordering is ex¬ 
pected to occur. 

Prospectively, it is of interest to incorporate interactions 
which will require more sophisticated theory but shall not 
change the basic mechanism of the measurement scheme. 
The superlattice modulation applied to a homogeneous 
system injects non-zero quasimomentum k^ into the sys¬ 
tem. Quite generally, one may realize arbitrary quasi- 
momentum transfer K if choosing different geometry of 
the perturbation, i.e. replacing the dimerization (— l) 3 
in Eq. (Al) by cos (Kj). This is promising as a spec¬ 
troscopic probe as it allows for the investigation of more 
complex lattice models and various intriguing quantum 
phases. 


ACKNOWLEDGMENTS 


We thank S. Wolff, A. Rosch and the group of M. Kohl, 
in particular L. Miller, J. Drewes and F. Brennecke, for 
insightful discussions. We acknowledge DFG (BCGS) for 
financial support. 


Appendix A: The superlattice modulation of the 
dimerized equilibrium lattice 


We diagonalize the dimerized tight-binding Hamilto¬ 
nian given by Eq. <§• First, we assign two sub-lattices 
with annihilation (creation) operators on odd sites and 
even sites defined by 


,(t) _ At) 

'jo — Wf+l.ffi 


&<!> = 


At) 

'2 j,a- 


Note that the first Brillouin zone is divided in half ] — 
fci/2, At/ 2] compared to the homogeneous equilibrium 
lattice. We then employ the Fourier transform, 


E< 


,ik2aj 


ttk.cn 


^E 7 E e ‘ fc2 % 


where k = k x , combined with a Bogoliubov transforma¬ 
tion, 

(3k,a = Ukdk,i7 T 


where 


u k = ^ ex P 


— I- 


,ip(k) — ka 


Vk = yQ exP 


. ip(k) — ka 


with 
ip(k) = 


I arctan 


(J' — J) sin(fca) 

( J'-\-J) cos (ka) 

b — k L 
Av 2 




With this we obtain the diagonal equilibrium Hamilto¬ 
nian Eq. (10) of the superlattice which exhibits two 
energy bands ±Ek = ± \J J’ 2 + J 2 + 2 J J' cos(2 ka). We 
approximate the superlattice perturbation by the dimer¬ 
ization operator 


o = E<-F( 

3,<? 


C j,cr C j+l,o A h- C ’ ) ■ 


(Al) 


Using the above transformation we obtain 

0 = J2 C °( fc ) (“Irak,* - Pt, <,&,*) 

k,cr 

+ C(k) (4,^,* - h.c. 


with 

Co(fc) = 2sin(fca) sin[^p(fc)], 

C(k) = 2isin(/co) cos[<p(fc)]. 

The resulting atom excitation rate to the upper band is 
given by Eq. where 


= co 2 

(j +J') 2 ] 

(j + j') 2 (j 2 + j' 2 ) 

16J 2 J' 2 

2J 2 J' 2 

1 

r (j + j') 2 {j 2 - j ,2 ) 2 i 


-2 

j 2 j i2 

5 


Is'MI = 


4 T 7/ A [u; 2 —4(J 2 +J' 2 )] 2 
4JJ V i 64 J 2 J' 2 


;/2 
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